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ABSTRACT 

We investigate the high-frequency radio spectra of twenty low-luminosity active 
galactic nuclei (LLAGNs) with compact radio cores. Our millimetre survey with 
the Nobeyama Millimetre Array (NMA) and analyses of submillimetre archival 
data that had been obtained with the Submillimetre Common User Bolometer 
Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT) reveal the fol- 
lowing properties. At least half of the LLAGNs show inverted spectra between 15 
and 96 GHz; we use published data at 15 GHz with the Very Large Array (VLA) 
in a 0.15-arcsec resolution and our measurements at 96 GHz with the NMA in 
a 7-arcsec resolution. The inverted spectra are not artificially made due to their 
unmatched beam sizes, because of little diffuse contamination from dust, Hii re- 
gions, or extended jets in these LLAGNs. Such high-frequency inverted spectra 
are apparently consistent with a 'submillimetre bump', which is predicted by an 
advection-dominated accretion flow (ADAF) model. We find a strong correlation 
between the high-frequency spectral index and low-frequency core power mea- 
sured with very-long-baseline-interferometry (VLBI) instruments. The inverted 
spectra were found exclusively in low-core-power sources, while steep spectra 
were in high-core-power ones with prominent pc-scale jets. This suggests that 
the ADAF and nonthermal jets may coexists. The flux ratios between disc and 
jet seem to be different from LLAGN to LLAGN; disc components can be seen 
in nuclear radio spectra only if the jets are faint. 
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1 INTRODUCTION 

Low- luminosity active galactic nuclei (LLAGNs), opera- 
tionally defined as a Ha luminosity of < 10*" ergs s~^, 
reside in about 40 per cent of nearby bright galax- 
ies, and occupy majority of the population of AGNs 
fao. Filippenko. fc Sargai^ Il997allbh . Its low luminosi- 
ties are thought to be caused by very low accretion rate 
on to super-massive black holes. Such accretion has of- 
ten been explained by the model of optically-thin disc 
such as an advection-dominated acc r etion flow (ADAF; 
| lchimarulll977l: iNaravan fc Yil 11994: lAbramowicz et alJ 
|l993), rather than optically-thick disc ('standard disc'; 
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IShakura fc Sunvaevll973) . It is predicted that the ADAF 
has a 'submillimetre bump' in its spectral-energy distri- 
bution (SED), while the standard disc is responsible for 
a big-blue bump in the SED of luminous AGNs such as 
quasars. Since a high brightness temperature of ~ 10^'^ K 
in radio bands is expected from the ADAF model, the 
first imaging investigations for accretion discs must be 
promising with future very-long-baseline-interferometry 
(VLBI) instruments that will provide micro-arcsec reso- 
lutions (^ ^rab^^jSh^^fara^^^^y urghvl lioo^ : lUlvestaJ 
I2OOOI : iMivoshi. Kameno. fc Falckd l2003l) Therefore, the 
observational evidence of the ADAF is crucial in the 
LLAGN radio sources. 

Although the ADAF model successfully explains 
the broadband spectra of LLAGNs, there is a signif- 
icant gap between observations and the ADAF spec- 
trum especially at low-frequency radio bands (e.g.. 
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[^faiim£to j_j>iineshige. fc Kusunos3ll997| : lOuataert et al.l 
This indicates that additional unknown 
emission does exist, and putative emission from the 
accretion disc may be buried under it. The submil- 
limetre bump means a highly inverted spectrum at 
centimetre-to-submillimetre wavelengths of spectral in- 
dex Q ~ .4, where Si, oc i^" *"", S,^ is flux density at fre- 
quency 1/ l|Mahadevanlll997t ). Observations for LLAGNs 
have been carried out exclusively at centimetre bands 
where high sensitivities are available, because most of 
LLAGNs are very faint radio sources. About half of 
low-luminosity Seyfert galaxies an d low-ionization nu- 
clear emission-line regions (LINER s;'Heckman 1980') hold 
a compact radio core l)Nagar et al. 2000. 2002a), and 
at leEist 25 per cent of transition objects, which are 
hypot hesized to be LINER/H ii-nucleu s composit e sys- 
tems JHo. Filinnenkn fc Sargend Eoo!^: lH;;] ri99^. also 
hold a compact core jFilho. Bart hcl. fc Ho 2000, 2002), 
at 15 or 8.4 GHz in 0.15-2.5 arcscc resolutions. Most 
of the cores show nearly flat spectra (a ~ 0) in 0.5- 
arcsec re solution: that is evidence for jet domination 
ijNagar. W ilson, fc Falcke 2001). Ahhough slightly in- 
verted spectra have been found in unresolved compact 
cores of several LLAGNs in milli-arcsec resolutions of 
VLBIs, their radio-to-X-ray luminosity ratios suggest 
that another significant radio component seems to con- 
tribute in_addition_toorigm the ADAF 
model (ni]vesta.d fc Hnl2001 ^lAnderson. TIlve:stad. fc Hri 
l2004) . To explain this rad io excess, an iet-ADAF 
model has been proposed ([Falcke fc BiermannI Il999l : 
lYuan. Markoff. fc Falck3 l2002l) . Thus, at centimetre 
bands, the contamination that is believed to be from jets 
has prevented the ADAF from being revealed. 

In the present paper, we report millimetre survey 
and submillimetre data analyses for many LLAGNs. Al- 
though technically difficult, a high-frequency observation 
is very promising to detect the inverted spectrum of the 
ADAF because of following two advantages: (1) spec- 
tral measurements at high frequencies are less affected by 
non-thermal jets, which generally show power-law spec- 
tra (-Q = 0.75-0.5), and (2) the ADAF predicts larger 
flux densities at higher frequencies. In fact, it has been 
obsevationally confirmed that flux densities at millime- 
tre wavelengths are 5-10 times large r than those at cen - 
timetre wavelengths in Sgr A* (e.g. iFalcke et alJll998j) . 
which is the nearest low-accretion rate system from us. 
However, contamination from diffuse dust may be harm- 
ful at > 100 GHz when we use a poor spatial resolution. 
Therefore, the use of a beam as small as possible and 
the estimations of every conceivable contamination are 
essential. 

The present paper is structured as follows. Sample 
selection is described in Section|5| Observations and data 
reduction are described in Section |3] In Section 01 we 
report these results and the estimations of diffuse con- 
tamination in our millimetre measurements. Our results 
of submillimetre photometry are utilized only for esti- 
mation for dust contribution. The origins of the spectra 
are discussed in Section |3 from the relation between the 
high-frequency radio spectrum and low-frequency radio 
core power. Finally, we summarize in Section |S| 



2 SAMPLE 

Our sample, 'VLBI-detected LLAGN sample', covers 20 
out of 25 all-known LLAGN radio sources that have 
been detected with VLBIs (column [1] of Table 0. 
The other five LLAGNs are N GC 4235, NGC 4450 
l|Anderson. Ulvestad. fc Ho 2004Y NGC 524, NGC 5354, 
and NGC 5846 jFilho et al. 20041) . which were recently 
newly VLBI-detected, but whose reports had not yet 
been published when we plan the survey. It is difficult to 
estimate resultant selection bias including into our sam- 
ple, because of multiplicity in selection criteria across the 
past VLBI surveys. However, at least we can say that 16 
out of 20 sources of our sample are all known LLAGNs 
that are at < 20 Mpc with > 2.7 mjy at 15 GHz, and the 
other four sources (NGC 266, NGC 3147, NGC 3226, and 
NGC 4772) are more distant. It may be, therefore, as a 
moderately good distance-limited and radio-fiux-limited 
sample for LLAGNs with a compact radio core. 



3 OBSERVATIONS AND DATA 
REDUCTION 

3.1 Millimetre band 

We made millimetre continuum observations for 16 out of 
20 targets using the Nobeyama Millimetre Array (NMA) , 
at the Nobeyama Radio Observatory (NRO), with D con- 
figuration that is the most compact array configuration. 
We had precluded the other four sources from our obser- 
vation list because they had already been observed at ~ 
3 mm in th e past (NGC 266, J)oi,_K ameno, & In ouc 2 005,; 
NGC 3031. ISakamoto et al.l2o'o¥NGC 4258, Doi et al. in 
prep. ; NGC 4486. iDespring^^^raix-Burnet. fc Davousli 
Il996l) . Our campaign spent more than 20 observational 
days between 2002 November 28 and 2003 May 25. Most 
of weak targets were observed over several days for inte- 
gration. Visibility data were obtained with double-sided- 
band receiving system at centre frequencies of 89.725 and 
101.725 GHz, which were Doppler-track ed. We used the 
Ultra Wide-Band Correlator (UWBC; lOkumura et alj 
I2OOOI) . which can process a bandwidth of 1 GHz per 
each sideband, i.e., 2 GHz in total. The wide bandwidth 
gives the NMA a very high sensitivity for continuum 
observations. Even if the systemic velocities of our tar- 
gets and galaxies' rotations of several hundred km s^^ 
are accounted for, the observing band can avoid possi- 
ble contaminations from several significant line emissions: 
"CO(J = 1-0), &^0{J =1-0), HCN(J =1-0), 
HCO+(J = 1-0), and SiO(J = 2 - 1). A system 
noise temperature, Tsys, was typically about 150 K in 
the double-sided bands. For observations of each target, 
we made scans of a reference calibrator close to the tar- 
get every 20 or 25 minutes for gain calibration. Bandpass 
calibrators of very bright quasars were scanned once a 
day 

The data were reduced using the UVP ROC-II pack- 
age iTsutsumi. Morita. fc Umevamalll997^ . developed at 
the NRO, by standard manners, including flagging bad 
data, baseline correction, opacity correction, bandpass 
calibration and gain calibration. To achieve higher sensi- 
tivity, visibilities of both the sidebands were combined in 
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the same weight, which resulted in a centre frequency of 
95.725 GHz. Each daily image was individually made in 
natural weighting and deconvolved using the AIPS soft- 
ware, developed at the National Radio Astronomy Obser- 
vatory (NRAO). Statistically significant variability was 
detected from the daily images in several LLAGNs, which 
will be reported in a separate paper. In the present pa- 
per, the visibilities of all days were combined by weighting 
with Tgys , and then imaged. The half-power beam widths 
(HPBWs) of synthesized beams were typically ~ 7 arcsec, 
corresponding to ~660 pc at a mean distance of 18.9 Mpc 
in our sample. The flux scales of the gain calibrators 
were derived with the uncertainty to 10 per cent by rel- 
ative comparisons to the flux-known calibrators, such as 
Uranus, Neptune, or Mars. These were scanned quasi- 
simultaneously when they were at almost the same ele- 
vations. The rms of noise on the images were estimated 
from statistics on off-source blank sky with the IMEAN 
task of the AIPS. Source identifications and measure- 
ments of flux densities were done on image domain by 
elliptical Gaussian profile fitting with the JMFIT task of 
the AIPS. We derived total errors in the flux measure- 
ments from root sum square of the errors in the Gaussian 
fitting (including thermal noise) and the flux scaling. 



in flux measurements from root sum square of the errors 
of the image noise and the flux scale. 



4 RESULTS 

Results of the NMA observations and the analyses of 
JCMT archival data are listed in Table Q We adopt a 
detection limit of 3cr""'' to peak intensities at the nu- 
clei. In detected cases, signal-to-noise ratios, S/N, are 
4.0-27 and 7.5-30 in the NMA and JCMT measure- 
ments, respectively. Although a spurious emission of ~ 3- 
4a™° sometimes appears in field of view of the NMA, 
we believe our detections because all of detected emis- 
sion components are found as a single point source at 
just the phase-tracking centre to an accuracy of less than 
HPBW/(S'/A^). We set the phase-tracking centres where 
compact components had been found in the past VLBI- 
astrometric observations. Because of S/N > 7.5 there is 
also no doubt in the detections with the JCMT. In unde- 
tected cases, the upper limit of flux densities takes 3o"™°. 

In addition to the past observations in literatures, we 
finally obtain centimetre-submillimetre continuum radio 
spectra for the LLAGNs of our sample (Fig. and their 
references are listed in Table HJ. 



3.2 Submillimetre band 

We attempted to measure flux densities at submil- 
limetre band from the nuclear regions using archival 
data. We searched the data that had been obtained 
at 347.38 GHz (850 fim) using Submillimetre Comm on 
User Bolometer Array (SCUBA; iHolland et 'aI]|l999^ on 
the James Clerk Maxwell Telescope (JCMT). The data 
of 14 out of 20 objects of the VLBI-detected sample 
were available to the public on 2004. The JCMT mea- 
surements for five objects have already been reported 
(NGC 266. lDoi. Kameno. fc Ino'u3l2005t NGC 4258, Doi 
etal. in prep.; N GC 4374. iLeeuw. Sansom. fc RobsonI 
| 2000l: NG C 4472, di Matteo et al.l Il999l: NGC 4486. 
iKnapp fc Patten 1991) . 

Four objects had been observed in jiggle-map mode. 
We reduced these data using the SCUBA User Reduc- 
tion Facility (SURF) package. Standard reduction pro- 
cedures were used, including flat fielding, flagging of tran- 
sient spikes, correction for extinction, pointing correc- 
tion, sky removal and flux-density scaling. Flux calibra- 
tors (Uranus, Mars, or CRL 618) gave us flux scaling 
factors (we assumed its uncertainty to ~15 per cent) and 
an effective beam size of 15.1 arcsec. The derivations of 
the rms of noise on images, source identifications, mea- 
surements of flux densities, and total error estimations 
were done in the same manners as the NMA reductions 
described in Section [3. II 

The other five objects had been observed in photom- 
etry mode. We followed standard reduction procedures, 
including flat fielding, flagging of transient spikes, cor- 
rection for extinction, sky removal, averaging, and flux- 
density scaling, using the SURF package. The photom- 
etry can measure a flux density from a beam-size region 
at the nucleus, i.e., an intensity. We derived total errors 



4.1 Radio morphology 

In the NMA observations, all of the detected sources are 
found as a single radio component. No significant ex- 
tended emission was found in the images with the NMA 
beam of ~ 7 arcsec. This indicates that there is no de- 
tectable emission of extended component in outer regions 
of > 500 pc from the nuclei. However, some unidenti- 
fied components other than the compact core could con- 
tribute into the NMA beam. This effect will be evaluated 
in Section 

On the other hand, in the JCMT measurements 
at 347 GHz, we found significant elongated features in 
NGC 4258 (Doi et al. in preparation), NGC 4579, and 
NGC 5866, and a very diffuse structure in NGC 3147 
(Fig. |5J . These features are consistent with the morphol- 
ogy of CO emissions that had been imaged with mil- 
limetric inte rferometers ( NGC 4258, iHelfer et al1l2003l : 
NGC 4579, ISofue et al.l l2003t) . with dust lane along 
the major axis of t he host galaxy (NGC 5866, e.g., 
i Malhotra et al.l l2000h . or with morphology of a low - 
frequency radio image (NGC 3147. lHummel et alJll985^ . 
It is, therefore, plausible that most of the 347-GHz emis- 
sion from these sources originates in stellar processes of 
star-forming regions, probably dust emission, associated 
with their host galaxies. 

4.2 Spectral indices of a 
centimetre-to-millimetre band 

Two-point spectral indices between 15 GHz (2 cm) and 
96 GHz (3 mm), which had never been investigated, are 
derived in this section. Unmatched beam sizes at differ- 
ent frequencies leads to incorrect spectral-index measure- 
ments, if significant diffuse components are present (a 
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Figure 1. Radio spectra of the VLBI-detected LLAGN sample. The symbols connected with solid lines represent quasi- 
simultaneous observations. For undetected data, 3ct™° upper limits are illustrated with downward arrows, where cr™'' is 
r.m.s. of image noise. Only for JCMT data at 347 GHz, intensities at nuclei, rather than total flux densities, are plotted. 
Dashed lines represent putative spectra only for visualizing the upper limits of potential contributions to measured 96-GHz 
flux densities from diffuse nonthermal (oc u~^'^), free-free (oc u~^'^), and dust (ex , where v is frequency) components 
(sec Section l4.4l and Tablelsl . Only for NGC 4143, no putative dust spectrum is not described because no data for estimation 
of dust contribution was available. Free-free emissions in NGC 4374, NGC 4472, NGC 4486, and NGC 4552 are too weak 
to appear in the ranges of the plots. 
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Table 1. Observation results for the VLBI-detected LLAGN sample. 



Name SgsGHz '^leGH. _ ■f347GHz_ 5347GHz <^3?7GHz_ ^ode Ref. 

(mjy) (inJy beam"-"^) (mJy beam~^) (mJy) (mJy bcam"-"^) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 266 


< 2.7 


0.9 


a 


< 18.5 


< 18.5 


6.2 


M 


a 


NGC 2787 


15.4 ±2.4 


1.0 




34.3 ± 5.9 




3.0 


P 


b 


NGC 3031 


390 ± 39*^ 




c 


86.6 ± 13.2 




2.1 


P 


b 


NGC 3147 


3.4 ± 1.5 


0.9 


b 


50.7 ± 9.3 


1747 ± 178 


5.4 


M 


b 


NGC 3169 


8.1 ± 4.0 


1.7 


b 












NGC 3226 


8.4 ± 1.8 


1.0 


b 












NGC 3718 


11.1 ± 3.5 


1.9 


b 


19.1 ± 3.8 




2.6 


P 


b 


NGC 4143 


6.7 ±2.7 


1.2 


b 












NGC 4168 


6.0 ± 2.5 


1.3 


b 












NGC 4203 


15.3 ± 2.9 


1.0 


b 


< 10.24 




3.4 


P 


b 


NGC 4258 


11.0 ±2.7*2 




d 


93.7 ±29.3 


3377 ± 254 


10.5 


M 


d 


NGC 4278 


57.0 ±8.5 


3.1 


b 


49.4 ±7.4 


53.3 ±8.0 


4.4 


M 


b 


NGC 4374 


147 ± 17 


5.0 


b 


144 ± 13*'* 






PM 




NGC 4472 


15.6 ±5.4 


2.1 


b 


3.0 ± 1.3 






P 




NGC 4486 


1880 ± 100*3 




e 


1069 ± 39 






P 




NGC 4552 


12.2 ±4.1 


2.6 


b 


< 26.6 


< 26.6 


8.9 


P 




NGC 4565 


< 4.7 


1.6 


b 












NGC 4579 


28.6 ±8.2 


2.7 


b 


85.7 ±23.8 


150 ± 30 


7.5 


M 


b 


NGC 4772 


7.1 ±2.4 


1.1 


b 












NGC 5866 


< 2.9 


1.0 


b 


51.7 ± 19.1 


173 ±31 


3.7 


M 


b 



*1 — mean flux density of multi-epoch data at 103 GHz. 

*2 — mean flux density of multi-epoch data at 101 GHz. 

*3 — obtained with the IRAM Plateau de Bure interferometer at 89 GHz. 

*4 — mean flux density of multi-epoch data. 



Note. — Column are: (1) galaxy name; (2) flux density at 96 GHz; (3) r.m.s. of image noise at 96 GHz; (4) reference 
for flux density at 96 GHz; (5) peak intensity of nuclear component at 347 GHz; (6) integrated flux density at 347 GHz; 
(7) r.m.s. of image noise at 347 GHz; (8) observation mode of the SCUBA at 347 GHz, 'M' represents jiggle mapping, 'P' 
represents photometry; (9) reference for data at 347 GHz. 



Note. — Reference: lallDoi. Kameno. & Inouf 


12003 b: the oresent DaDer:ld:ISakamoto et al.l200ll; d: Doi et al. in DreD.:la: 


iDesorinffre. Fraix-Burnet. &; Davoustll99fiHi 


Leeuw. Sansom. &; R.obson 




:ldi Matteo et al.ll99fll:lhHKnaDD Hr, PattenI 


Il99ll 



'resolution effect'). Extended jets are one of the candi- 
dates. Because of their power-law spectra this sort of con- 
tamination would be small at millimetre, but possible at 
centimetre wavelengths if we use a large beam. Hence, as 
the lower-frequency data t o derive the two-poi i it spec tral 
indices, we use data from iNaear et al.l ll200d. l2002all of 
the Very Large Array (VLA) A-configuration at 15 GHz 
with a resolution of ~ 0.15 arcsec, corresponding to 
~ 14 pc at 18.9 Mpc, which is a mean distance of our 
sample. These are unique and uniform data set that can 
provide centimetre flux densities at the highest frequency 
and in the highest spatial resolution for all our LLAGN 
sample. The other contaminations such as dust emission 
and free-free emission from H ll regions could contribute 
especially to the NMA beam. Very careful estimations 
for these diffuse contributions will be described in Sec- 
tion lOl 

The two-point radio spectra are listed in column (6) 
and (7) in Table |21 These results are from unoperated, 
measured flux densities at 15 and 96 GHz, for which 
the diffuse contaminations are not removed. We adopt 
a threshold of a = between inverted and steep spectra. 
The threshold definitely distinguishes 13 out of 20 sources 
beyond the margins of error of their spectral indices. The 
other 6 out of remaining 7 sources were marginally dis- 



tinguished; the threshold lies within the margins of er- 
ror of their spectral indices. We classified them into four 
spectral categories: definitely inverted T, marginally in- 
verted 'i', marginally steep 's', and definitely steep 'S' 
spectra. Consequently, we found inverted spectra in many 
LLAGNs. However, a substantial fraction of the LLAGNs 
showed steep spectra. 

4.3 Submillimetre emission 

Submillimetre spectra of compact cores are a very inter- 
esting issue for observational examination of the ADAF 
model. However, the JCMT fiux densities could be domi- 
nated by dust emission from host galaxies, since beam 
widths of the JCMT were too large. If the core was 
brighter than the dust, it may be detectable from its 
flux variation even with the large beam. NGC 4374, 
NGC 4472, NGC 4579, and NGC 5866 had been observed 
at multiple epochs with the JCMT. A signiflcant vari- 
ability was seen in NGC 4579: 85.7 ± 23.8 mJy beam"^ 
on 2002 January 16 by jiggle-map mode and 25.8 ± 
4.9 mJy beam~^ on 2001 March 15 by photometry mode. 
The significant fraction of observed submillimetre fluxes 
may be direct emission of the core. In the other sources, 
no significant variation was found. At least, the results 
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Figure 1 - continued 



from the JCMT measurements should be useful for es- 
timation of the upper limits of dust contribution (Sec- 
tion mUTl 



4.4 Diffuse fractions in 96-GHz flux densities 

The beam mismatch between the measurements of 15 and 
96 GHz might pollute the two-point spectral indices. Be- 
cause the NMA beam at 96 GHz, ~ 7 arcsec, was greatly 
larger than the VLA one at 15 GHz, ~ 0.15 arcsec, the 
derived spectra would be biased toward being inverted 
if some diffuse components exist in addition to the core. 
Possible contaminations are (1) dust emission, (2) free- 
free emission from plasma heated/ionized by UV radia- 
tion from OB stars in nuclear star-forming regions or an 
accretion disc of LLAGN, and (3) extended non-thermal 
emission from jets or supernova remnants. In the next 
three subsections, we estimate the contaminations from 
the three diffuse components to the NMA measurements. 



4-4-^ Contamination from dust emission 

The dust emission at a given frequency can be estimated 
by extrapolating an appropriate spectrum that was deter- 
mined by measurements with the same aperture at higher 
frequencies where the dust contribution predominates 
any other components. The graybody spectrum of dust 
emission is S^"'^'^ oc u^^'^ in the Rayleigh- Jeans regime 
at less than about 1 THz, where (3 is the dust omissiv- 
ity, w hich is believed to be between 1 and 2 jHildcbrandl 
From the JCMT resuhs at 347 GHz, we can de- 
rive useful upper limits of the dust contribution to the 
observed 96-GHz flux densities by following three overes- 
timates. 

(i) The spectral indices of the dust spectra differ from 
galaxy to galaxy. The l3 value has been determined to 
be 1.3 ± 0.2 by JCMT observations for nearby galaxies 
jPunne et al.. 2000'). To keep it safe, we adopt (3 = 1, 
which is the most influential case in the spectral extrap- 
olation to lower frequencies. This must provide an over- 
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NGC3147 




DEC 




10h17m06s 10h17m00s 10h16m54s 10h16m48s 10h16m42s 
RA 



12h20m08s 12h20m06s 
RA 





RA 



15h06m33s 15h06m30s 
RA 



Figure 2. JCMT images at 347 GHz for all the VLBI-detected LLAGNs that had been observed in SCUBA jiggle-map mode, 
excep t for NGC 266 (Doi. Kamcno. & Inouc 200^), NGC 4258 (Doi et al. in prep.), and NGG 4374 Sansom, R obsoiJ 



NGC 5866, 



"347GHz ' ' ^""^ -t^iu, i^347Qjjz' 

T™^^Qjj^ x3,4,5,. . . and 15, where o-™7^qj]2 r.m.s. of image noise. The HPBWs of the JCMT, 15.1 arcsec, are 
illustrated at the lower left-hand corner of the images. Astromctric core positions are described in cross symbols. NGC 3147, 
NGC 4579, and NGC 5866 are spatially resolved. 



estimate for dust contribution at 96 GHz. (ii) Beam size 
of the JCMT, 15.1 arcsec, is about two times larger than 
that of the NMA, so it must lead to also an overesti- 
mate as photometry for the NMA aperture, (iii) Contri- 
butions from the other components such as a core, jets, 
Hll regions, and molecular clouds [CO(J = 3 — 2) lines at 
345.796 GHz] etc., could contribute the JCMT measure- 
ments, so it must also make an overestimate compared to 
a measurement only for dust emission. 

Thus, from the three overestimates, we can derive 



an upper limit of the contamination from dust in the 
NMA measurement for each LLAGN. For sources with- 
out the JCMT measurement, we used flux-density data 
from literatures: a single-dish observati on at 230 GHz 
with a 20-arcsec beam for NGC 4565 (|Neininge^^^ 
Il996r) : the IRAS Faint Source Catalogue llMoshi^^d 
|l99Cl) at 60 and 100 fim with ~l-2 arcmin beams and 
the relati on between far-infr ared (FIR) and 850-/im to- 
tal fluxes iPunne et alj|2000h for NGC 3169, NGC 3226, 
NGC 4168, and NGC 4772. Unfortunately, no appropri- 
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Table 3. Source parameters for the VLBI-detected LLAGN sample. 



Name 


Class 


Dist. 
(Mpc) 


Ref. 


Morph. 


15GHz 
96GHz 


Sd 


Ref. 


cdust 

f mJv) 


'-'max 


cnonth 
'-'max 

('mjv'l 


rnonth 
J max 


j-total 
J max 


(1) 




(3) 


(4) 


(5) 


(6) 


(7) 


f8) 


(9) 


(10) 








NGC 266 


L1.9 


62.4 


a 


SB(rs)ab 


< -0.23 


s 




0.4 


0.2 


1.0 






NGC 2787 


LI. 9 


7.5 




SB(r)0+ 


0.42±0.14 


I 




0.7 


0.7 


1.3 


0.09 


0.18 


NGC 3031 


SI. 5 


3.6 


c 


SA(s)ab 


0.46±0.09 


I 




1.8 


9.7 


10.4 


0.03 


0.06 


NGC 3147 


S2 


40.9 


z 


SA(rs)bc 


-0.47±0.50 


s 




1.1 


0.6 


10.9 


3.18 


3.68 


NGC 3169 


L2 


19.7 


T 


SA(s)a pec 


0.09±0.57 


i 




14.6*1 


1.0 


10.5 


1.30 


3.22 


NGC 3226 


LI. 9 


23.6 


T 


E2: pec 


0.24±0.21 


I 




11.1*1 


0.5 


0.4 


0.05 


1.43 


NGC 3718 


LI. 9 


17.0 


a 


SB(s)a pec 


0.01±0.34 


i 




0.4 


0.4 


1.8 


0.16 


0.23 


NGC 4143 


LI. 9 


15.9 


T 


SAB(s)0 


0.39±0.41 


i 






0.6 


1.2 


0.18 




NGC 4168 


SI. 9 


30.9 


T 


E2 


0.35±0.43 


i 




0.5*1 


0.0 


0.7 


0.11 


0.21 


NGC 4203 


LI. 9 


15.1 




SABO-: 


0.26±0.19 


I 




0.2 


0.8 


0.7 


0.05 


0.11 


NGC 4258 


S1.9 


7.3 




SAB(s)bc 


0.70±0.21 


I 




2.0 


2.2 


0.7 


0.07 


0.44 


NGC 4278 


LI. 9 


16.1 




E1+ 


-0.25±0.15 


S 




1.0 


5.5 


46.6 


0.82 


0.93 


NGC 4374 


L2 


18.4 




El 


-0.12±0.11 


S 




3.0 


1.0 


785*^ 


5.36 


5.38 


IN LjV^ 44 ( Z 


GO 
OZ 


1 ^^ Q 


u 

L> 




U. ( ZitU.OD 


T 
1 




L). i 


U.U 


Oft ft 


i. ( U 


1.(1 


NGC 4486 


L2 


16.1 




EO + pec 


-0.22±0.03 


S 




22.4 


3.4 


16748 


8.91 


8.92 


NGC 4552 


T2 


15.3 




E 


-0.85±0.38 


S 




0.6 


0.4 


12.1 


0.99 


1.07 


NGC 4565 


SI. 9 


17.5 




SA(s)b?spin 


<0.13 






1.4*2 


0.7 


6.7 






NGC 4579 


S/L1.9 


19.1 


e 


SAB{rs)b 


0.01±0.32 


i 




1.8 


3.7 


11.7 


0.41 


0.60 


NGC 4772 


LI. 9 


16.3 


a 


SA(s)a 


0.39±0.34 


I 




0.5*1 


0.4 


0.1 


0.02 


0.15 


NGC 5866 


T2 


15.3 




SAO + spin 


< -0.52 


S 




1.1 


0.1 


2.6 







*1 — derived from IRAS data, see text in detail. 
*2 — derived from 230-GHz data, see text in detail. 

*3 — derived from data obtained with a single-dish telescope, see text in detail. 

Note. — Column are: (1) gala xy name; (2) nuclear activity type as determined in Palomar optical spectroscopic survey 
Filippenko. fc Sargentlll997a^ . 'L' represents LINER, 'S' represents Seyfert, 'T' represents objects with transitional L + Hll 
region-type spectrum. '2' implies no broad Ha is detected, '1.9' implies broad Ho is present but not broad H/9, and '1.5' implies that 
broad Ha and broad H/3 a re detected; (3) distance to galaxy ; (4) reference of the distance, as listed below; (5) Morphology type 
of host galaxy, as listed in IhoT Filippenko. &: SargentI il997al) : (6) two-point spectral index between 15 GHz and 95 GHz, derived 
from data with the VLA A configuration and NMA, respectively; (7) spectral category, classified by a threshold of OggQ^^ = 0. 'S' 
and 's' represent steep spectra, ctggQ^^ < 0, determined definitely and marginally, respectively. T and 'i' represent inverted spectra, 
'^seGHz = ^' determined definitely and marginally, respectively; (8) reference for the VLA data at 15 GHz, as listed below; (9) upper 
limit of the potential contribution at 96 GHz from dust emission. This was estimated mainly from peak intensity at 347 GHz and by 
extrapolating a dust spectrum, see text in de tail; (10) upper limit of the potential contribution at 9 6 GHz from free— free emission. 
This was estimated from H/3 fiux measured bv lHo et alJ ll997tf) and lHo. Filippenko. fc fiargenj ^^2Z3)i see text in detail; (11) upper 
limit of the potential contribution at 96 GHz from diffuse nonthermal emission. This was estimated from data at 1.4 GHz, see text 
in detail; (12) upper limit of diffuse fraction only of the nonthermal component in the 96 GHz fiux density; (13) upper limit of total 
diffuse fraction in the 96 GHz flux density. 

Note. — Referenc e: 13 lTullvlll9'8l Ibt iTonrv et aljl200ll: S iFreedman et alJll994t |3 iHerrnstein et al.lll999t S ISolanes et ai]l2002l: Q 
iNagar et alJl200(i l3 lNagar et alJl2002al 



ate measurement for NGC 4143 at higher frequencies 
than 96 GHz was available. The upper limits of the flux 
densities of the diffuse dust emission at 96 GHz, SjJiaxi 
are listed in column (9) of Table 01 These values suggest 
that there must be little influence to the NMA measure- 
ments for all the LLAGNs, except for a few sources that 
was evaluated using the IRAS data. The IRAS beams 
may be too large only for the nuclear region. Therefore, 
all the NMA measurements are probably little affected 
from dust contamination. While we adopted the extreme 
of dust spectra with /3 = 1, there is nevertheless the slim 
possibilities that we may have detected at 96 GHz an un- 
usual dust componen t with an abundance of very cold or 
very small grains (cf. iLisenfeld et 'ai]l2002h . 



Contamination from free-free emission 

Free-free emission from Hll regions in nuclear star- 
formation or narrow- /broad-line regions around the AGN 
can be independently estimated from the extinction- 
corrected flux of H/3 line emission by assuming an 

electron temperature and metallici ty of the plasma 

JCaplan fc DeharvendlTgsl ICondonI Il993) . If the elec- 
tron temperature is ~ 10'' K and iV(He+)/iV(H+) ~ 0.08, 

(f)—i;;i^)(Gk) 

where is the expected flux density of free-free emis- 
sion at frequency u and F(H/3) is Hf3 line flux. The ex- 
tinction values and extinction-corrected fluxes of broad 
and narrow H/3 for our sample can be derived from 
the data of the Palomar optical spectroscopic survey 
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Table 2. Reference list for Fig. 1. 



Name 



(1) 



VLA-D 



(2) 



VLA-A/ 
VLA-B 

(3) 



VLBl NMA/ 
JCMT 

(4) (5) 



NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 



266 

2787 

3031 

3147 

3169 

3226 

3718 

4143 

4168 

4203 

4258 

4278 

4374 

4472 

4486 

4552 

4565 

4579 

4772 

5866 



ab 
a 

ab 
ab 
ab 

^ 

ab 

& 

ab 

ab 
a 



be 
dc 

eb 

edc 

ebc 

edc 

ebc 

edc 

ebc 

dc 

edc 

dc 

edc 

ebc 

ebc 

edc 

ebc 



T 

hi 

1 

ki 

dm 
nm 
T 

T 



xr 

X 
X 
X 



*1 — mean flux density of multi-epoch data. 
*2 — data obtained with the IRAM Plateau de Bure interfer- 
ometer at 89 GHz. 

Note. — Column arc: (1) galaxy name; (2) reference of data 
obtained with the VLA D configuration at 1.4 or 15 GHz; (3) 
reference of data obtained with the VLA A or B configuration 
at 1.4-15 GHz; (4) reference of data obtained with VLBI; (5) 
reference of data obtained with the NMA and JMCT. 

Reference: [a| [Condon et al 



Note. 



20001: 13 iNagar. Wilsom fc Falckell2001 
" Becker. White, fc Helfandl llgSd 
J l200d: H It 



ippgt Ibt iNagar et al.l 



Nagar et al.l2002, 



a. 



iBartel fc Bietenholz 



lUlvestad 



Falcke et al.l 120 
Hoi l2001al: 



mi!! 



K ameno. fc Inou g [ 200,4 IT lAndg 

IIP 



Hcrrnstcin et al 
■ 

Lv. Walker. &: Wrobo 



1997t II J 

J200I S iFomalont et al.l [ 2000 

3|Safa^^toIet_alJ |200lj; p: Do i et al. in prep.: 
Leeuw. Sansom. fc RobsonI EOOCI: H Idi Matteo et 



19991: 1st iDesprin^^^Sx-Burnet. fc Davousr " ^9i 
KnapD fc Pattenlll99]r x: the present paper 



||Ho. FiliDDenko. & Sareend 


1997al:lHoetal.lll997cD us- 


ing_the extinction curve of 


Cardelli. Clayton. & Mathid 



il989fl . We can derive useful upper limits of the free-free 
contribution to the NMA measurements by following two 
overestimates. 

(i) Amount of hidden broad line components have 
been unknown in the optical survey. In its sample, frac- 
tions of broad li ne flux in total Ha flux are at most 0.98 
\Ro et al.lll997cl) . which may be the least hidden and the 
most broad-line-dominant case. We adopted this frac- 
tion, then the upper limit of the broad line fluxes can 
be derived using extinction-corrected narrow line fluxes. 
It may be an overestimate of an intrinsic total line flux, 
(ii) This optical survey was done with 2x4 arcsec^ aper- 
ture, ~ 4.8 times smaller than that of the NMA 7-arcsec 
beam. Hence, we magnified the extinction-corrected nar- 
row line H/3 flux from a 2 x 4 arcsec^ aperture to 4.8 times, 
then obtained an equivalent flux to the 7-arcsec region. 



It can provide an overestimate because we assumed here 
that the surface density of narrow-line clouds is constant 
over the 7-arcsec region, although that tends to be rather 
highly central concentrated (J-'ogEC ct al. 200^. We dealt 
here with only narrow-line components, which could be 
distributed over about a few hundred parsecs from the 
nucleus, while broad-line ones should concentrate within 
a sub-parsec region. 

From the two overestimates, we can derive the upper 
limit of the free-free emission at 96 GHz, Smax, in the 
NMA measurement for each LLAGN. The results, listed 
in column (10) of Table It suggests that there is little 
influence to the NMA measurements. 



4.4-3 Contamination from extended non-thermal 
emission 

The extended non-thermal components of jets or super- 
nova remnants must be optically thin, so show steep spec- 
tra with Q ^ —0.5. A lower-frequency observation with 
a larger beam should be more sensitive to this emission. 
With published data that had been obtained with larger 
beams at 1.4 GHz, we can deduce useful upper limits 
of contamination from the extended non-thermal compo- 
nents to the NMA measurements by following the two 
overestimates. 

(i) We use data of th e NRAO VLA Sky Survey 
fNVSS: ICondon et al.ll99l at 1.4 GHz with ~ 45-arcsec 
beam, which is signiflcantly larger than that of the NMA. 
For two sources, NGC 4258 and NGC 4772, uncovered in 
the NVSS, we used data of the Faint Images of the Ra- 
dio Sky at Twenty-centimetre (FIRST) survey with ~ 5- 
arcsec beam. The beam of the FIRST survey was slightly 
smaller than that of the NMA, but emission with a size 
of ~ 7 arcsec was expected to be detectable without be- 
ing resolved out. Only NGC 4374 was u ncovered in both 
the catalogues, therefore we used data of lWhite fc Becker! 
{1992) with ~ 700-arcsec beam. These 1.4 GHz measure- 
ments give us estimations of amount of extended non- 
thermal components in roughly the same or larger region 
than the NMA beam. It makes an overestimate of a non- 
thermal flux into the NMA beam, (ii) We adopt a = —0.5 
that would give the most influential case by extrapolat- 
ing a power- law spectrum from 1.4 to 96 GHz. It makes 
an overestimate of a nonthermal flux density at 96 GHz. 

From the two overestimates, we can deduce an upper 
limit of the extended nonthermal contamination, S'JJia'x''i 
in the NMA measurement for each source. The results 
are listed in column (11) of Tabled For about half of 
our sample, very large contamination is possible. Most of 
them show jet-like structures, which have been found by 
past observations (Sections 15.11 and 15.211 . For the other 
half, little contamination is ensured. 



4-4-4 Total diffuse contaminations and spectral indices 

For our NMA measurements, we obtain the upper lim- 
its of possible contaminating flux densities at 96 GHz, 
SmTx, Smax, and S^ax"^, from dust, free-free, and non- 
thermal emission components, respectively. The upper 
limit of total contamination fraction is now defined as 
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ft. 

J IT 



[S^ax + + «'']/S96GHz. The derived val- 



ues are listed in column (13) of Table El The value of 
Smax**^ accounts for the major part of /max' for almost all 
the LLAGNs, except for a few sources whose dust con- 
tributions was estimated using the IRAS data with the 
much larger beam sizes. This can be seen by comparing 
/maS*'' With where /"Tx'^ = SZxV^QecHz listed 

in column (12) of Table|21 It seems that both the diffuse 
dust and free-free emission components are practically 
negligible in arcsec-resolution observations at ~ 3 mm 
for LLAGNs. It is only necessary to concern that the 
diffuse nonthermal components could be responsible for 
artificially-made inverted spectra between 15 and 96 GHz 
due to unmatched resolutions; the LLAGNs with large 
/max' would tcud to be biased toward inverted spectra. 

However, our results do indicate rather an inverse 
tendency. Almost all the LLAGN showing inverted spec- 
tra have a practically negligible diffuse contamination 
with /^i°x'' < 0.5, while almost all the LLAGNs show- 
ing steep spectra have possibly large contaminations with 
/m'S"' > 0.5 (Fig. OJ. A Wilcoxon test about the spec- 
tral indices between the LLAGNs with /max'" < 0.5 
and those with f^^^ > 0.5 resulted in a probability of 
0.0246 for the null hypothesis that these two groups were 
drawn from the same parent population. Note that three 
LLAGNs whose f^°^^ or aggQH^ values cannot had been 
be determined were excluded in this statistical test. 

Thus, for almost all the LLAGNs showing inverted 
spectra, we cannot find any evidence that the observed 
two-point spectra between 15 and 96 GHz had been ar- 
tificially made by resolution effect. The observed spectra 
presumably reflect the intrinsic spectra of the cores that 
had been detected as unresolved sources with the VLA 
in a 0.15-arcsec resolution. Also, possible flux variability 
between observational epochs at 15 and 96 GHz could 
not have artificially made such many inverted spectra un- 
less the spectra were intrinsically inverte d; the LLAGNs 
sometimes show rapid variability (e.g., iHo et al.l Il999l : 
ISakamoto et al.ll200ll:fNagar et al.ll2002a^ . but our result 
must be statistically supported by a large number of our 
sample. 



5 DISSCUSSION 

5.1 Origins of radio spectra 

We found inverted high-frequency radio spectra in a num- 
ber of LLAGNs with compact cores. On the other hand, 
a significant fraction of our sample showed steep spectra 
fSection l4.21 . We discuss the variety of the observed radio 
spectra. An interpretation is that optically-thin nonther- 
mal jets may be responsible for these steep spectra, while 
the inverted spectra may be caused by ADAF discs. The 
presence of jets are supposed to be recognized by very 
high radio powers that cannot be generated in the ac- 
cretion disc, or by its characteristic radio morphology. In 
this section, we investigate the radio core power and ra- 
dio morphology of our sample, in order to have clues to 
the origins of the observed radio spectra. 

The relation between the two-point spectral in- 
dex (agecHz) and the 5-GHz core power {PsaSI) that had 




I 

(a>o) 



s 

Spectral category 



Figure 3. Histogram of two-point spectral indices between 
15 and 96 GHz, ai^GHz^ f^j. ^j^g VLBI-detected LLAGN sam- 
ple. Filled region represents sources with f^^^ < 0.5, where 
= 5Si°'i*'V 596GHz (see Section HX!!! and and Ta- 
ble|2j. It is warranted that diffuse contamination is practically 
negligible for these sources. Open region represents sources 
with fup^^^ 0.5; we cannot ruled out significant diffuse con- 
tamination for these sources. Grey region represent sources 
whose diffuse fractions cannot be estimated because of unde- 
tection in their NMA measurements. NGC 4565 is not included 
here because of its deficient determination of spectral category. 
This histogram confirms that there are many LLAGNs show- 
ing inverted spectra that are not artificially made by resolution 
effect. 



been measured in ~ 1 milli-arcsec resolutions of VLBIs 
(column [2] in Table |1J is shown in Fig.0] The informa- 
tion about jet morphology (column [3] and [5] in TablellJ 
is also included in this plot. From this plot, we can find 
the following three properties: (1) a clear inverse correla- 
tion, (2) there is no LLAGN showing inverted spectrum 
in the region of high core power, and (3) the LLAGNs 
with pc-scale jets are found exclusively in the region of 
steep spectrum and higher core power. 

An inverse correlation between the spectral index 
and VLBI core power can be clearly seen as Fig. |1] The 
LLAGNs with higher VLBI-core powers tend to show 

steeper spectra. A generalized Kendall's r test, u 

the ASURV s tatistical package flsob e fc Feigelsonlll99i 
iLaVallev. Isob c. & Fcigolson 1992), said that a probabil- 
ity is only 0.0023 for the null hypothesis that there is 
no correlation between them. This suggests that the ob- 
served spectra are strongly related to the VLBI core pow- 
ers. The VLBI core powers are expected to include both 
the contributions from compact jets and the accretion 
disc. However, because of the relatively low frequency, 
5 GHz, the core power is expected to be biased toward the 
jet contribution. On the other hand, because of relatively 
high frequencies, the spectral index may be biased toward 
the disc contribution. The inverse correlation might ap- 
pear as a consequence of the combination of such biases. 
In other words, the inverse correlation suggests an intrin- 
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sic nature of LLAGNs: j ets and disc coexist in the cor e 
(a jet-ADAF model, e.g.. lYuan. Markoff. fc Falck J2002I ) : 
the observed steep spectra are as the result of jet domina- 
tion, and the observed inverted spectra are as the result 
of disc domination. 

There is no LLAGNs showing inverted spectra in 
the region of high core power that cannot be generated 
by th e ADAF. In an simple ADAF model ijMahadevanI 
Il997l) . the maximum radio power at a given frequency is 
mainly dependent on black h ole mass and accretion rate 
(cf. t he same discussion in lAnderson. Ulvestad. fc Hoi 
l2004r) . With a black hole mass of 10® Mq, the ADAF 
cannot generate a 5-GHz radio power greater than ~ 
8.9x10^" WHz"^ in any accretion rate (as a vertical line 
in Fig|lJ. It is consistent with the apparent upper limit of 
the VLBI core power of LLAGNs showing inverted spec- 
tra. Hence, we can reasonably say that the ADAF may 
be responsible for the inverted spectra for the lower-core- 
power objects, while potential disc components may be 
dominated by jet components for the higher-core-power 
objects. 

All the four LLAGNs with pc-scale jets are found ex- 
clusively in the region of steep spectrum and higher core 
power, as filled symbols in Fig. |1] Note that the 15-GHz 
flux densities, which were used for derivation of the two- 
point spectral indices, were from a region of <0.15 arc- 
sec, typically corresponding to < 14 pc from the nuclei, 
where the pc-scale jets have been found. It is, therefore, 
not surprising that the four LLAGNs were categorized 
as definitely-steep spectra 'S' (Table 2J- Their high core 
powers are presumably generated by the base of these pc- 
scale jets, rather than discs. On the other hand, six out 
of nine LLAGNs with kpc-scale jets reside in the range 
of inverted spectrum. For four out of the six, we stress 
that their inverted spectra were not artificially made by 
resolution effect due to the kpc-scale jets because we have 
confirmed that the diffuse nonthermal fiux densities into 
the NMA measurements were significantly small (Sec- 
tion ^^31- That is, the six LLAGNs have both compact 
cores showing inverted spectra and very faint kpc-scale 
jets. Our NMA beam size was small enough to separate 
off their kpc-scale jets. For the other two LLAGNs with 
kpc-scale jets, NGC 4472 and NGC 3169, we cannot rule 
out resolution effect. The two LLAGNs are seen as the 
two open regions on the spectral category T and 'i' in 
Fig- El respectively. 

Consequently, we propose an idea that the high- 
frequency radio spectra mainly depend on which is more 
dominant in the core, jets or disc. This idea provides 
a reasonable explanation on the clear inverse correla- 
tion between the high-frequency radio spectrum and low- 
frequency core power. The inverted spectra of the ADAF 
may be observable only when the jet components are rela- 
tively faint. We may have detected direct radio emission 
from the accretion disc of the LLAGNs with faint jets. 
The power ratios between disc and jet seem to differ from 
object to object. It is very interesting for astrophysics on 
the issue of the accretion budget allocation on a black 
hole system. 

However, we must note that our results do 
not rule out the possibility that the inverted spec- 
tra are made only by jets. In a jet model pro- 



Table 4. Monochromatic 5-GHz core power ranking and 
known jets in the VLBI-detected LLAGN sample. 
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(1) 



1 pVLBI 
^"S ^5GHz 
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(2) 



Jet 
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Ref. 
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Jet 
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Ref. 



(6) 



Sp. 
(7) 



4ZOo 


1 n 


Q 
O 


o aD 


Iv 


o^ 


T 

1 


IN CjO 4uz 


19 


6 








T 
1 


IN KjLj 44 / Z 


19 


9 


c 


IV 


PQ 


T 
1 


TVjpp 9787 


1 Q 


Q 








J 


INf-jL.- 40Di3 


on 


i 


a 








-\J/~i f~i Q'71 O 


20 


3 




iv 


r 


i 


iN<jO OODU 


20 


4 








c 

b 


NKjtLj 4zUv> 


20 


4 








T 
i 


NGC 3031 


20 


4 


S Kh 


K 


s 


I 


NGC 4143 


20 


4 


c 






i 


NGC 4168 


20 


4 


rf 






i 


NGC 3169 


20 


5 


z 


K 


t 


i 


NGC 3226 


20 


5 








I 


NGC 4579 


21 





dfe 


K 


a 


i 


NGC 4278 


21 


1 


P dji 


K 


u 


S 


NGC 3147 


21 


2 








s 


NGC 266 


21 


2 


dk 






S 


NGC 4552 


21 


4 


P ~ 






s 


NGC 4374 


21 


8 


P ci 


K 


V 


s 


NGC 4486 


22 


5 


P Imn 


K 


w 


s 



Note. — Column are: (1) galaxy name; (2) core power de- 
tected in VLBI observations at 5 GHz, except for NGC 4258 at 
22 GHz; (3) presence of jet feature in the VLBI image. 'S' rep- 
resents sub-pc-scale jets. 'P' represents pc-scale jets; (4) ref- 
erence for the VLBI observation, as listed below; (5) presence 
of jet feature in VLA image. 'K' represents kpc-scale jets; (6) 
reference for the VLA observation, as listed below; (7) spectral 
category, the same as Tablc|21 The sources showing steep spec- 
tra have exclusively high VLBI-core powers, > lO'^^ W Hz~^. 

Note . — Reference: lal 
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posed by iFalcke fc BiermannI il999|) . the submillime- 
tre bump is produced by the acceleration zone of 
the jet, called nozzle, associated with the ADAF, 
while the low-frequency fiat spectrum comes from the 
outer part of jet. The jet model can explain the ob- 
served high-luminosities of low-frequency radio emis- 
sion, which cannot be generated by a simple ADAF 
(lUlvestad fc Hoi l2001al: lYuan. Markoff. fc Falck^ I2OO2I : 
lAnderson. Ulvestad. fc Hdl2004h . However, the outer jet 
structures must have been revealed in the VLBI images 
if the jet axes were not fairly clo se to our line of sight 
(lAnderson. Ulvestad. fc Hdl2004ll . To judge it, imaging 
analyses in micro-arcsec resolutions would be necessary. 
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logPs^GHz [W Hz ^] 

Figure 4. Relation between 5-GHz radio power of core de- 
tected witli VLBIs and spectral index between 15 and 96 
GHz for the VLBI-detected LLAGN sample. Squares repre- 
sent sources with kpc-scale jets, while circles represent sources 
without kpc-scale jets. Filled ones represent sources with pc- 
scale jets f Table 111 . The dot-dashed line represents the radio 
power limit that the ADAF on a black hole with lO'' Mq can 
generate. 



5.2 Jet-disk symbioses in LLAGNs 

We have suggested in the previous section that jets 
and the ADAF seem to coexist in LLAGNs. However, 
only an unresolved structure has been found in milli- 
arcsec im ages for mos t of the V LBI-detecte d LLAGNs 
jFalcke et al. 2000: Ulvestad fc Ho 2Q01a; .Nagar et al.l 
l2002al: lAnderson. Ulvestad. fc H0II2OO4I 1. On the other 
hand, pc/sub-pc-scale jet structures have been found in 
several LLAGNs. We find from the core power ranking in 
Table|l]that the presence of VLBLdetectable pc-scale jets 
depends on the core radio powers. The pc-scale jets have 
been found exclusively in high-core-power sources. In low- 
core-power sources, no pc-scale jet has been detected in 
VLBI ima ges even with moderately high dy namic ranges 
of ~ 100 ilAnderson. Ulvestad. fc Hdl2004l . Cannot the 
low-core-power sources generate high-brightness jets? 

The cases of NGC 3031 and NGC 4258 are 
very suggestive. Sub-pc-scale jets h ave been re- 
vealed in their nuclei (e.g., JH crrnst ein et all Il997l : 
iBietenholz. Bartel. fc RupenI lioooll . nevertheless, they 
are low-core-power LLAGNs. Note that the two LLAGNs 
are the nearest two sources in the VLBI-detected sample 
(Table 121. We speculate that all low-core-power sources 
also have sub-pc-scale or smaller jets, and that such struc- 
tures would be revealed if they were at less than sev- 
eral Mpc like NGC 3031 and NGC 4258. The detectabil- 
ity of small-scale jets might merely depend on distance. If 
that is true, the 5-GHz VLBI-core power should contain 
jet contributions significantly. 

Additionally, NGC 3031 and NGC 4258 are very 
suggestive cases also about coexistence with the ADAF 



disc. Although they show clear evidence of jets, their 
high-frequency spectra are highly inverted, aggcHz = 
0.46 ±0.09 and a^^gg^ = 0.70 ±0.21, respectively. These 
are the highest two spectral indices in the VLBI-detected 
LLAGNs with negligible diffuse contaminations. If their 
millimetre emissions originate in the ADAF discs, the 
two LLAGNs would be examples for the coexistence of 
the jets and ADAF. 

We mention apparent exceptions. NGC 266 and 
NGC 3147 show steep spectra, but have no detectable 
jet in spite of their high core powers (Table |1J . Since 
the two LLAGNs are the most two distant LLAGNs in 
our sample, it could be somewhat hard to reveal their 
putative small-scale jet structures. Their high radio pow- 
ers might be also responsible for the dominance of such 
small-scale jets. NGC 5866 show a definitely steep spec- 
trum, in spite of a low core power. Its nuclear activity 
type, 'transition', might be related to the jet domina- 
tion; another transition object, NGC 4552, also shows a 
definitely steep spectrum. However, our limited sample of 
transition objects cannot allow us to discuss it in detail. 



5.3 Conflict with the unified scheme? 

There appears a clear difference in the spectral properties 
between type Is (SI. 5, SI. 9, and LI. 9) and type 2s (S2, 
L2, and T2), as seen in Tabled Most of the type Is show 
inverted spectra, while most of type 2s show steep spec- 
tra. We suggest this diff erence conflicts with the unifled 
scheme dAntonucc i^'iggS*) . In this scheme, many observed 
properties are determined basically only by a viewing an- 
gle with respect to the obscuring dusty torus surrounding 
the nucleus. Type-1 Seyferts, with broad permitted opti- 
cal emission lines, are those in which the torus is viewed 
closer to face-on orientation. Type-2 Seyferts lack broad 
permitted lines, except in polarized flux, and thought to 
be seen with the line of sight passing through the torus 
in an orientation closer to edge-on. This classification has 
been expediently used in LINER s and transitions (e.g., 
iHo. Filippenko. fc Sargen^ ll997b'l. Since the radio emis- 
sion is transparent through the dusty torus, luminosities 
and spectral indices of type Is and 2s should look just 
the same to the observer unless affected by relativistic 
beaming. However, radio powers of type 2s tend to be 
significantly larger than those of type Is in our sample. 
It seems that the difference is due to larger contribution 
of jets in type 2s, nevertheless, which would be less af- 
fected by relativistic beaming than type Is. 

There are also the other reports a bout differences 
between type Is and type 2s in LLAGNs. lUlvestad fc Hoi 
J200ld) found that low-luminosity type-1 Seyferts rather 
show somewhat stronger centi metre radi o emis sion than 
low-luminosity type-2 Seyferts. iHo et af] (120011) have re- 
ported that type-2 LLAGNs are underluminous in an X- 
ray band compared to type-1 ones with the same Ha 
luminosity. There are 'intrinsically true type-2 Seyferts 
(without broad lines even i n polarized fiux)' espec ially 
in the low-luminosity range fTran' 2003; Laor 2003. and 
reference therein). Thus, the discrepancy in the unified 
scheme is commonly seen in LLAGNs. 

We suggest a possible idea for the discrepancy 
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with the unified scheme in terms of radio properties of 
LLAGNs. If outfiows carry accreting mass away and de- 
crease a density at the innermost disc, then the spectral- 
peak frequency could decrease to several tens of GHz 
(|di Matteo et al]l200(t ). Although such a central engine 
could generate strong jets, neither high-frequency in- 
verted spectra nor a lot of up-scattered photons that 
will ionize broad-line clouds can be produced. That is, 
LLAGNs with strong jets would tend to have weak broad- 
line emission, which could not have been detected in a 
limited sensitivity of the optical survey, then they could 
be classified into type 2s, even if the obscuring torus is 
in a face-on orientation. 



coexist in the core of LLAGNs, and their flux ratios 
are different from object to object. It suggests that 
relative radiative efficiencies between the jet and disc 
would determined by some kind of physical parame- 
ter(s) in the central engine. From the theoretical point 
of view l|Nara^M^gunMnshch ev. fc Abramowic3 I2OO0I : 
iBlandfon^^^egehnemFQQOh . it has been sxiggested that 
radiatively inefficient fiows of low accretion rate tend 
to stimulate strong convective instabilities and power- 
ful outfiows. The central engines of LLAGNs probably 
can produce jets. We will discuss the physical parame- 
ters that control jet productions in our LLAGN sample 
in a separate paper (Doi et al. in prep). 



5.4 Implication for the ADAF 

We have found evidence for inverted high-frequency ra- 
dio spectra in a number of LLAGNs. Such a high- 
frequency property is potentially the long-sought sig- 
natures of the ADAF. It is an apparent observational 
confirmation for the submillimetre bump theoretically- 
predicted in the ADAF model. In past observations, 
the radio spectra of LLAGNs had been investigated 
exclusively at cent imetre bands where hig h sens itivi- 
ties are available. iNa gar. Wilson, fc Falckel (|200ll ) re- 
ported the radio spectra of 16 LLAGNs that had been 
observe d with the VLA at 5-15 GHz; lUlvestad fc Hoi 
reported the radio spectra of three LLAGNs 
that had been observed wit h the V LBA at 1.7-8.4 GHz; 
lAnderson. Ulvestad. fc Hoi (1200411 reported the radio 
spectra measured with the VLBA for three LLAGNs at 
2.3-15 GHz and three LLAGNs at 1.7-43 GHz. The fre- 
quency range of submillimetre-bump hunting has greatly 
expanded by our high-frequency radio survey. Our find- 
ings will encourage further LLAGN surveys at higher- 
frequency ranges with, e.g., the coming Atacama Large 
Millimetre and submillimetre Array (ALMA), to prove 
the ADAF. 

The existence of the ADAF will give radio as- 
tronomers a chance to make the direct imaging for ac- 
cretion discs. Their area of study had been limited to 
the nonthermal phenomena in relativistic jets of AGNs, 
although optical and X-ray astronomers have observed 
the direct fiuxes from the accretion discs as unresolved 
sources. Radio astronomers have VLBI techniques, which 
provide the greatest spatial resolutions all over the elec- 
tromagnetic wavelengths. A high-frequency VLBI instru- 
ment will report the first image of morphology of ac- 
cretion discs in micro-arcsec resolutions. The VSOP- 
2, a next generation space VLBI mission currently be- 
ing planned, will provide a 38 micro-arcsec resolution 
at 43 GHz and a brightness-temperature sen sitivity of 
10* K llHirabavashi. Murata. fc Murphvll2002ll . Because 
an electron temperature of io^~^-^ K at the region within 
several tens Schwarzschild radius is predicted by the 
ADAF model, the VSOP-2 can reveal disc structures of 
more than a dozen nearby LLAGNs. If the accretion disc 
were an optically-thick standard disc, it will have never 
been detected with the VLBI because of low brightness 
temperature. 

From our results, both the jets and disc seem to 



6 SUMMARY 

Our millimetre survey and analyses of submillimetre 
archival data for the 20 VLBI-detected LLAGNs reveal 
high-frequency radio spectra of LLAGNs with compact 
cores as follows: 

(i) At least half of the sources show inverted spectra 
between 15 and 96 GHz. We used the published data at 
15 GHz with the VLA in a 0.15-arcsec resolution and our 
measurements at 96 GHz with the NMA in a 7-arcsec 
resolution. 

(ii) Submillimetre fiuxes measured with the SCUBA on 
the JCMT may be from the dust associated with host 
galaxies. These data have been useful to estimate dust 
contamination in the NMA measurements. 

(iii) We confirm little diffuse contamination in the NMA 
measurements from dust, free-free, and extended non- 
thermal emissions for LLAGNs showing inverted spectra. 
Therefore, the inverted spectra are not caused by resolu- 
tion effect, but intrinsic properties of the cores. 

(iv) In our sample of 20 LLAGNs, radio cores with in- 
verted high frequency radio spectra reside in at least ten 
LLAGNs. 

(v) We find an inverse correlation between the high- 
frequency radio spectrum and the low-frequency radio 
core power. Such a property can be explained by an idea 
that fiux ratios between the jet showing steep spectrum 
and the ADAF disc showing inverted spectrum are differ- 
ent from LLAGN to LLAGN. The disc components could 
be seen only if the jets are faint. 

Consequently, we have found evidence for inverted high- 
frequency radio spectra in a number of LLAGNs. Such 
a high-frequency property is potentially the long-sought 
signatures of the ADAF. 
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